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motion  throughout  the  atmosphere. 

Case  studies  are  presented,  using  an  enhanced  version  of  the  Barnes  (1985) 
computer  program.  The  first  case  study  deals  with  a  severe  weather  situation  that 
resulted  in  tornadic  activity  through  eastern  Texas.  The  second  case  study  covers  a 
two  day  period  in  which  a  Mesoscale  Convective  Complex  (MCC)  moved  through  the 
central  plains  states.  Results  of  these  case  studies  show  that  the  Q-vector  form  of  the 
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the  initialization  of  the  prognostic  numerical  models.  I  he  diagnostic  charts,  from 
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however,  the  diagnostic  results  are  best  used  when  incorporated  with  other  routinely 
available  charts.  A  drawback  to  the  program  is  that  noiCgeostrophic  effects  such 
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IX  1  K'ODI  (TIO\ 

Before  ill**  operational  meteorologist  can  make  In-  next  ion-cast  lie  must 
examine  all  the  information  available  ami  then  benin  tin-  process  of  det.-nmnimz 
'.'.'h;it  factors  will  lx-  responsible  for  the  weather  durtne  the  forecast  period.  S« <n i«* 
■  «f  'lie  information  available  comes  in  the  form  of  computer  ueueraled  fun-casts  Mr 
different  at mospheric  processes.  Oik  of  the  important  prognostic  products  is  tin- 
vertical  velocit  v  chart  P  r«  xl  need  1  >  V  numerical  models  sill'll  as  the  Nested  (irid  Model 
f  N  ( ;  M  )  and  its  predecessor,  t  he  Limited  area  l-ine  mesh  Model  l  IT. Mi. 

\ertical  motion  in  tlx-  atmosphere  has  a  pronounced  effect  upon  the  'seatin'!'. 
I'pward  motion  produces  clouds,  cyclones,  fronts,  precipitation,  and.  in  general, 
destabilization  of  tin*  atmosphere.  Downward  motion  produces  anticyclones  and. 
in  general,  clearing  conditions,  flie  current  prognostic  proiluct  depicting  vertical 
motion  is  tin*  TOD  mb  vertical  velocity  chart.  A  question  here  is,  what  can  be  done 
to  improve  the  methods  for  evaluating  atmospheric  vertical  motion. 

I  he  classical  omega  equation  is  one  possible  diagnostic  equation  that  can  be 
used  to  determine  vertical  motion.  The  two  forcing  terms  for  vertical  motion,  in 
the  omega  equation,  are  1  I  the  vertical  derivative  of  vorticiiy  advection  <uul  2i  the 
horizontal  l.aplacian  of  temperature  advection.  Altlioiiuh  these  two  terms  can  be 
associated  with  seperate  physical  processes,  they  often  ha\e  a  degree  of  cancellation 

I  he  journal  style  used  is  from  Journal  of  tfn  I  tmoipht  nr  Sen  ma  s. 
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hrtuvf'n  tluMii.  Pin*  t « >  a  srrirs  of  approximation?-  in  tL<  vainm  «>l  the  uiim^a 
equation,  tin1  result  ant.  vertical  motion  is  applicable  <ud\  when  I  lie  atmosphere  jv, 
in  a  state  of  hvdrostatic  and  geostrophic  balance.  I  lie  lesnlts  are  best  used  in  the 
middle  t  r<  >posphere. 

In  the  l!)  Ill’s  the  evaluation  of  vertical  motion  equations  was  approached  troru 
two  different  viewpoints  according  to  Hoskins  et  ah  il9,Si.  I  hey  stated  that 
Sutcliffe  produced  his  tlieorv  of  (level: pment  from  a  diagnostic  point  of  view  while 
C’harney  extemh'd  his  tlieorv  of  harochnic  instabilitv  id\uaiuic  p<iint  of  novi. 

Sutcliffe  (  1917),  in  his  development,  neglected  terms  involving  the  deformation 
of  t  lie  wind  field.  These  terms,  as  Sut  cli  ffe  point  ed  out .  can  be  important  in  front  a  I 
regions  as  well  as  in  jet  entrance  and  exit  regions.  I  herelore.  in  large  deformation 
regions  such  as  frontal  zones  Sutcliffe's  approximations  are  le-s  applicable. 

As  stated  above,  the  classical  omega  equation  litis  a  drawback  in  that  tin' 
vorticity  and  temperature  terms  can  be  misleading  when  isolated  in  order  to 
determine  their  relative  magnitude  or  the  resulting  sign  of  vertical  motion.  1  his 
is  partly  due  to  the  vorticity  term  being  a  vertical  derivative  and  the  temperature 
term  being  determined  on  a  horizontal  plant*.  The  two  terms  are  not  independent. 
Trenbert  h  (  1 97S)  showed  that  the  vorticity  and  tempera)  n re  terms  can  he  expressed 
with  a  common  t  erm  between  them.  I  his  common  term  is  the  ad  vert  ion  of  vort  i  <i  t  y 
bv  tin'  thermal  wind.  In  the  vorticity  term,  the  common  term  is  a  positive  qmmtitv 
while  in  the  temperature  term  it  is  a  negative  quantity.  I  Inis,  the  two  terms  have 
a  degree  of  cancellation  between  them. 


In  this  lli*  ■sis  a  modified  version  of  the  omega  diagnostic  equation  wd!  be 
presented  along  the  lines  indicated  bv  Hoskins  and  his  collaborators  (Hoskins. 
Draghici.  and  Davies  1978)  that  overcomes  these  problem  areas.  It  will  provide 


a  determination  of  vertical  motion  fields  based  upon  the  divergence  of  a  quantity 
called  the  Q-vector,  A  three  level  numerical  model  (Barnes  1985)  that  uses  this 
new  vertical  motion  equation  will  be  discussed.  Case  studies  using  the  model  will 
be  shown  in  order  to  determine  the  model’s  usefulness  as  an  enhanced  tool  for 
operational  meteorologists. 

The  usefulness  will  be  partly  determined  by  making  two  sets  of  comparisons. 
The  first  comparison  will  be  to  see  if  precipitation,  as  depicted  by  the  National 
Weather  Service  (NW7S)  Radar  Summary  chart,  is  occurring  in  a  area  of  computed 
upward  vertical  motion.  This  will  be  done  by  overlaying  the  Radar  Summary  charts 
onto  the  800,  600,  and  400  mb  V  •  Q  charts.  A  second  comparison  will  be  to  see 
how  much  of  the  computed  subsidence  areas  is  precipitation  free.  Both  sets  of 
comparisons  will  be  overlaid  onto  a  grid  marked  in  1°  increments.  Results  in  each 
set  of  comparisons  will  be  shown  as  a  set  of  percentages  for  each  level. 

Durran  and  Snellman  (1987)  as  well  as  Doswell  (1987)  have  presented  argu¬ 
ments  for  the  use  of  the  Q-vector  form  of  the  omega  equation.  These  authors  agree 
that  it  is  the  best  way  to  evaluate  vertical  motion  numerically. 
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CHAPTER  II 


VERTICAL  MOTION  EQUATIONS 


Hoskins,  Draghici,  and  Davies  (1978,  hereafter  HDD)  published  a  new  deriva¬ 
tion  of  the  omega  equation  that  included  the  terms  that  Sutcliffe  omitted  in  his 
derivation.  They  derived  a  one-term  representation  for  the  geostrophic  determi¬ 
nation  of  ageostrophic  motion.  From  their  derivation,  the  synoptic  distribution  of 
vertical  motion  can  be  determined  from  a  height  and  temperature  field  as  given 
on  a  constant  pressure  chart.  Hoskins  and  Pedder  (1980)  demonstrated  that  this 
new  method  of  deriving  vertical  motion  can  provide  more  information  than  the  Sut¬ 
cliffe  equation  in  describing  details  of  synoptic  system  development  especially  when 
dealing  with  the  horizontal  ageostrophic  field  and  intense  frontal  circulations. 

The  derivation  given  by  HDD  used  a  vertical  coordinate  z  —  (B^Cp/  g)(l  - 
(p/Po)*)-  The  basic  equations  are: 

du  d4> 

~dt  ~  V  ~  dx  ’ 

dv  d(f) 

It  = 

de  n 
di~ °’ 

du  dv  dw 
dx  dy  dz 

“d  r,  =  (i)6 

where  60  and  Po  are  standard  values  of  potential  temperature  and  pressure,  R  is  the 
gas  constant,  «  =  R/CP,  f  is  the  Coriolis  parameter  and  assumed  to  be  constant, 


KW 


and  is  the  geopotential.  All  other  symbols  represent  their  standard  meteorological 
meanings. 

The  geostrophic  wind  components  are: 

1  d<t>  ,  1  d<f> 

u9~~~rz~  and  v9-~rz~- 

f  dy  f  dx 

The  thermal  wind  relations  are: 

dug  _  g  d6  dvg  _  g  d$ 

f  dz  ~  e0  dy  and  f  dz  ~  60dx' 

The  geostrophic  vertical  component  of  relative  vorticity  is: 

=  (2) 

Taking  the  vertical  derivative  of  (2),  and  then  substituting  the  geostrophic 
terms  from  (1)  into  the  result  yields 


—  iL  |  I  \ 

1  dz  e0  \dx*  dy2)  e0 

Now  using  the  quasi-geostrophic  form  of  the  vorticity  equation, 


A  +  v  ,vV  =/— 
at  9  r9  1  dz' 


taking  its  vertical  derivative,  multiplying  through  by  /,  and  substituting  from  (3) 


produces  the  equation 


The  quasi-geostrophic  and  adiabatic  form  of  the  First  Law  of  Thermodynamics 


is  given  as 


i-  +  vs  •  6  =  -w-r© 

at  9  J  dz 


(5) 


,  ■%.  *  >*  »  ■  * 


where  @(z )  is  the  potential  temperature. 
Multiplying  (5)  by  (g/60)V2  yields 


|-|-V20  +  |-V2(Vfl  •  V0)  +  N2V2w  =  0 

6q  at  Gq 


where  N 2  =  ( g/6o)(d,Q/dz )  is  the  squared  buoyancy  frequency  and  is  assumed  to 
be  independent  of  x  and  y  in  this  derivation. 

Eliminating  the  time  derivative  between  (4)  and  (6)  produces 

A'^w  +  =  /A(V,  ■  VC,)  -  |-V’(V,  •  V«),  (7) 

which  is  the  usual  form  of  the  omega  equation.  This  equation  shows  that  the  vertical 
motion  is  a  result  of  the  vertical  derivative  of  vorticity  advection  and  the  horizontal 
Laplacian  of  thermal  advection. 

In  his  derivation  of  a  vertical  motion  equation,  Sutcliffe  neglected  adiabatic 
cooling  and  warming;  therefore  the  N2V2w  term  did  not  appear  on  the  left  hand 
side  of  (7).  The  forcing  on  the  right  hand  side  of  the  equation  was  reduced  to  one 
term  at  the  expense  of  part  of  the  horizontal  Laplacian  of  thermal  advection. 

A  new  version  of  the  omega  equation  is  derived,  as  described  in  HDD,  by  first 
multipling  the  second  equation  of  motion  through  by  /  and  taking  the  vertical 
derivative.  The  result  is 

9  (  fdVA  I  V  Vf9v»  ,  fdV9  V7r  ,  .  n 

dt{fJ^  +y’  vfJI  +  ,~aT'Vv’  +  f  nr  ~  °" 


Now  the  term 


=  fdUg  dVg  ,  f9Vg  dVg 

dz  9  dzdx  dz  dy 


can  be  rewritten  using  the  continuity  equation  as 


,dVg  dug  dv  dvgdu 

f  dz  '  V  ~  f  dz  dx  ~  f  dz  dx' 


By  substituting  the  thermal  wind  components  (1)  it  follows  that 


,8V  g  r- 

f-dT'Vv 


Therefore  the  second  equation  of  motion  can  be  written  as 


i+v..v)/£  +  <?1+^-0. 


(8) 


Using  the  thermal  equation  and  multiplying  through  by  g/60  and  then  taking 
the  x  derivative  the  result  is 


d  (9  de\  ,  g  8V g  ^  ,  g -T  -86  g  86  dw 


J..X .  v—  |  y  =  0 

8t  \60  8x  J  60  8x  60  9  8x  6 0  8z  8x 


The  following  substitutions  can  be  made: 


9_  86~  _  ,8vg 
60  8x  8z 


V6  =  -Qj, 


9  dVg 

60  8x 


g  86  dw 
60  8z  8x 


and 


dw 

8x 


Thus  the  thermal  wind  equation  can  be  rewritten  as 


(9) 


In  (8)  it  is  seen  that  the  vertical  shear  of  the  ageostrophic  wind  component, 


uag,  makes  a  contribution  to  local  change  of  the  horizontal  gradient  of  potential 
temperature.  In  (9),  the  horizontal  gradient  of  potential  temperature  is  affected  by 
the  horizontal  variation  of  w. 
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Eliminating  the  time  derivative  in  (8)  and  (9)  by  subtracting  (8)  from  (9)  yields 

S**  ■  (10) 

As  seen  in  (10),  Qi  depends  on  the  horizontal  variation  of  vertical  motion  and 
the  vertical  shear  of  ageostrophic  motion  in  the  x  direction.  This  helps  to  maintain 
the  thermal  wind  balance. 

By  similar  arguments  it  can  be  shown  that 

2Q,  =  (11) 

Oy  Oz 

Talcing  the  x  derivative  of  (10),  the  y  derivative  of  (11),  and  adding  the  two 
resultant  equations  produces  the  vertical  velocity  equation 


where 


+  =  2V-Q 


The  left  hand  side  of  (12)  is  identical  to  the  left  hand  side  of  (7)  which  is  the 
usual  form  of  the  omega  equation.  The  right  hand  side  is  now  a  single  forcing  term 
commonly  referred  to  as  the  divergence  of  the  Q-vector. 

Quasi-geostrophic  vertical  motion  is  a  result  of  the  tendency  for  advection  by 
the  geostrophic  w'ind  to  destroy  thermal  wind  balance.  In  reality  the  atmosphere  is 
never  really  far  from  thermal  wind  balance.  The  balance  is  continually  reestablished 
by  ageostrophic  circulations.  Hoskins  and  Pedder  (1980)  showed  that  the  Q-vector 
provides  an  approximate  picture  of  the  ageostrophic  horizontal  wind  in  the  lower 
branch  of  the  circulation  that  develops  in  order  to  maintain  thermal  wind  balance  in 
evolving  synoptic  disturbances.  In  quasi-geostrophic  theory,  on  an  f-plane  (Coriolis 
parameter  is  held  constant)  the  vertical  velocity  is  determined  by  the  divergence  of 


„  4-.  I.-  v  r  ^  /  «-  r.  V.  S  «-  «•  /.  -V  V  -"  V  f 
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PHYSICAL  INTERPRETATION  OF  Q 


In  order  to  see  the  physical  significance  of  Q  we  can  assume  a  rectangular 
cartesian  coordinate  system  with  the  j  axis  tangential  t.»  the  potential  temperature 
at  the  point  where  the  value  of  Q  is  to  he  determined.  'The  ;/  axis  will  point  towards 
the  colder  air.  In  this  rase  8d  i). r  -  0  and  the  (.)- vector  components  hecome 


and  (}>  — 


9  ‘h'y  89 
dtl  ()j-  ()y 
g  di'q  i)d 
do  Ou  <hj 


No w  is  related  to  the  horizontal  shear  and  (),  to  the  confluence  or  difHuence 
of  the  geostrophic  motion.  This  is  illustrated  in  Fig.  1  and  Fig.  2.  As  was  stated 
above,  the  potential  temperature  is  only  a  function  of  y  in  both  of  these  figures. 


F I f 5 .  1.  Q,  horizontal  shear  case.  Here  v-  i.-,  increasing  with  x.  Ihe  heavy 
arrows  represent  Q  and  show  that  is  positive  while  (),  (I.  1  he  y  component 

of  motion  (thin  arrows)  is  positive  at  A  and  negative  at  IF 


I  he  magnitude  of  Q  in  each  case  is  inversely  proportional  to  the  spacing  of  tin* 
potential  temperature  contours  and  directly  proportional  to  the  gradient  of  the  v 


FIG.  2.  Q,  confluent  case.  In  this  case  vg  is  decreasing  with  y.  The  heavy 
arrows  represent  Q  and  show  Qi  =  0  and  that  Q 2  is  negative.  The  y  component  of 
motion  (thin  axrows)  is  negative  in  the  cold  air  and  positive  in  the  warm  air. 

component  of  geostrophic  motion.  This  can  be  written  as  Q  a  (l/L)X?vg  where  L 
is  the  spacing  of  the  6  contours. 

From  the  original  omega  equation  u>  >  0  implies  descent  while  u>  <  0  implies 
ascent.  Equation  (12)  then  shows  that  V-Q  >  0  corresponds  to  descending  motion 
while  V-Q  <  0  is  related  to  ascending  motion.  Assuming  that  Q  is  negligible 
outside  points  A  and  B  then  in  Fig.  1,  at  point  A,  the  northward  moving  air  should 
be  ascending  while  the  southward  moving  air,  at  point  B,  is  descending.  If  we  were 
to  draw  streamlines  on  the  figure  then  we  would  see  Q  in  a  trough.  Thus  we  have 
rising  motion  ahead  of  the  trough  and  subsidence  behind  the  trough.  As  is  the 
case  in  the  omega  equation,  rising  motion  can  be  associated  with  positive  vorticity 
advection  and  subsidence  with  negative  vorticity  advection. 

Figure  2  shows  a  confluent  case  with  vg  decreasing  with  y.  In  this  figure  it  is 
assumed  that  Q  is  negligible  to  the  north  and  south  of  the  Q-vectors  drawn.  Then 
this  figure  shows  that  the  southward  moving  air  is  subsiding  in  the  cold  region  and 
is  associated  with  divergence  of  Q  while  the  northward  moving  air  is  ascending  in 
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the  warm  sector  (V  •  Q  <  0). 

In  each  case  the  Q-vector  is  pointing  toward  the  area  with  rising  motion.  This 
is  also  shown  in  Figs.  3  and  4  which  depict  frontogenetic  and  frontolytic  situations 
respectively.  Figures  3  and  4  assume  that  the  magnitude  of  the  Q-vector  is  negligible 
to  the  left  and  right  sides  of  the  Q-vector  drawn. 


WARM 


FIG.  3.  Frontogenetic  circulation.  Slanted  lines  represent  a  frontal  boundary. 
Heavy  arrows  represent  Q.  Small  arrows  represent  the  ageostrophic  circulation. 


WARM 


FIG.  4.  Frontolytic  circulation.  Slanted  lines  represent  a  fronted  boundary. 
Heavy  arrows  represent  Q.  Small  arrows  represent  the  ageostrophic  circulation. 

In  Fig.  3  the  geostrophic  field  tends  to  intensify  the  temperature  gradient  and 
the  direction  of  Q  is  toward  warmer  air.  The  low  level  ageostrophic  flow  that  is 
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established  to  maintain  the  thermal  wind  balance  is  towards  the  warm  side.  The 
cross  frontal  circulation,  in  this  case,  has  the  warm  air  rising  ahead  of  the  front 
and  cold  air  subsiding  behind  the  front.  This  is  a  typical  frontogenetic  circulation 
pattern. 

Figure  4  shows  the  opposite  case  with  the  temperature  gradient  being  weakened 
and  the  cross  frontal  circulation  implying  frontolysis.  There  is  warm  air  subsiding 
and  cold  air  rising. 

It  may  be  noticed  that  in  both  Fig.  3  and  Fig.  4  that  the  Q  is  pointing  in 
the  same  direction  as  the  low  level  ageostrophic  motion.  Having  a  chart  of  Q  then 
gives  a  vectorial  representation  of  horizontal  ageostrophic  motion  below  the  level  at 
which  Q  is  computed.  The  same  chart  then  gives  some  indication  of  the  intensity  of 
the  fronted  circulation.  Having  a  map  of  Q  superimposed  over  a  temperature  field 
can  then  show  where  V6  is  being  enhanced  or  weakened  and  thus  give  a  picture  of 
system  development. 

From  the  above  discussion  it  can  be  seen  that  having  a  chart  depicting  Q  (or 
a  chart  of  V  •  Q)  shows  more  than  just  vertical  motion.  We  can  more  readily  see 
frontal  circulations  and  their  corresponding  effects  as  to  whether  a  front  has  the 
tendency  to  strengthen  or  not.  The  Q-vector  includes  terms  that  were  neglected 
by  vertical  velocity  charts  currently  in  use.  Thus  the  Q-vector  form  of  the  omega 
equation  gives  a  more  complete  representation  of  the  quasi -geostrophic  rorcing  of 
vertical  motion  than  the  traditional  isobaric  maps. 
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A  THREE  LEVEL  Q-VECTOR  PROGRAM 


Cahir,  Norman,  and  Lowry  (1981)  applied  the  Q- vector  formulation  to  a 
summer  time  cold  front  moving  across  the  Great  Lakes.  This  was  done  as  part  of  a 
demonstration  of  a  real  time  computer  graphics  sj'stem  that  could  be  implemented 
for  operational  forecasting  purposes.  The  demonstration  only  implemented  the  Q- 
vector  formulation  for  the  700  mb  level  and  the  authors  found  that  the  ascent  forced 
by  quasi-geostrophic  dynamics  was  spatially  related  to  the  cold  front. 

Barnes  (1985)  carried  the  Q-vector  formulation  further  by  introducing  a  two- 
level  model.  He  used  reported  height  data  at  mandatory  pressure  levels  to  compute 
contoured  patterns  of  vertical  motion  and  other  diagnostic  tools  at  the  two  levels. 
The  model  was  applied  to  a  severe  convective  storm  over  eastern  Colorado  that 
resulted  in  the  development  of  a  mesoscale  convective  complex  (MCC). 

For  the  Barnes’  case  the  LFM/MOS  and  quantitative  precipitation  guidance 
indicated  that  eastern  Wyoming  and  western  Nebraska  as  the  most  threatened  area, 
whereas  the  actual  convection  was  most  intense  in  Colorado.  The  resulting  MCC 
propagated  into  Kansas  and  weakened.  Another  convective  system  developed  in 
northeastern  New  Mexico  and  moved  into  the  Texas  panhandle  while  intensifying. 
Barnes’  Q-vector  diagnostics  correctly  showed  the  favorable,  for  storm  development, 
large  scale  forcing  over  eastern  Colorado  and  also  over  northeastern  New  Mexico  4 
h  before  the  observed  changes  in  the  convective  systems.  According  to  Barnes,  the 
LFM  guidance  showed  that  the  Kansas  storm  would  have  the  most  favorable  large 
scale  forcing  for  development  and  thus  was  in  contradiction  with  his  results. 
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Barnes  has  since  enhanced  his  model  to  produce  three  levels  of  contoured 
diagnostic  maps.  Through  personal  communication  a  micro-computer  version  of 
his  model  has  been  obtained.  The  program  takes  approximately  8  min,  on  a  8MHz 
IBM  AT  class  machine,  to  produce  20  diagnostic  maps.  This  time  does  not  include 
the  data  ingest  time  which,  for  this  study,  takes  about  5  min  using  a  modem  to 
download  data  from  the  Texas  A&M  Meteorology  Department’s  Harris  computer. 
The  required  data  for  the  model  are  the  850,  700,  500,  and  300  mb  geopotential 
height,  temperature,  and  dew  point  depression. 

The  remainder  of  this  section  will  deed  with  the  structure  of  Barnes’  three  level 


model  and  some  of  the  significant  algorithms  used.  The  results  of  using  this  model 
will  be  presented  in  the  case  study  section  of  this  thesis. 

The  Barnes’  program  produces  much  more  than  maps  of  Q,  V  •  Q,  and 
a  geostrophic  frontogenetic  function  (GFF).  GFF  is  the  tendency  (moving  with 
the  geostrophic  motion)  for  frontogenesis  or  frontolysis  to  occur  and  is  evaluated 
by  determining  the  component  of  the  Q-vector  along  the  potential  temperature 
gradient  (i.e.,  GFF  oc  Q  ■  V0).  The  model  also  produces  lapse  rate  of  temperature, 
mixing  ratio,  height  contours,  temperature  advection,  thickness,  vorticity,  and 
different  types  of  stability  maps  among  others.  At  the  user’s  discretion,  20 
additional  maps  can  be  contoured  and  the  program  can  easily  be  modified  to  produce 
other  desired  maps. 

The  program,  as  written  by  Barnes,  is  in  BASIC  and  can  be  rewritten  in 
FORTRAN  if  the  need  arises.  The  only  problem  may  be  in  rewriting  the  graphics 
output.  Linkable  graphics  routines,  such  as  the  National  Center  for  Atmospheric 
Research  (NCAR)  graphics  package,  would  have  to  be  used  since  FORTRAN  doesn’t 
implement  any  primitive  graphics  routines.  The  description  of  the  model  is  partly 
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extracted  from  the  information  received  with  the  program  and  also  from  ” dissecting” 
the  provided  code. 


The  model  was  developed  to  perform  calculations  on  a  portion  of  the  National 
Meteorological  Center’s  (NMC)  LFM  grid  that  is  used  at  the  Techniques  Develop¬ 
ment  Laboratory  (TDL).  The  LFM  grid  is  cartesian  with  AX  =  AY  =  190.5  km 
oriented  along  105°W  longitude.  The  origin  of  the  grid  is  near  the  southwest  corner 
of  the  North  American  continent  (see  point  1,1  in  Fig.  5).  The  TDL  grid  is  a  31  by 
24  portion  of  the  LFM  grid  and  begins  at  LFM  point  (15,  8).  The  model  uses  a  18 
by  15  portion  of  the  grid  and  the  position  of  the  grid  is  user  selectable.  The  grid 
and  map  background  must  be  selected  before  running  the  model  for  the  first  time 
or  if  a  new  grid  portion  is  desired.  The  origin  of  the  diagnostic  grid  is  the  northwest 
corner  and  the  north-south  grid  points  are  numbered  decreasing  northward. 

The  map  background  is  polar  stereographic.  For  output  to  the  screen  the 
coordinates  of  political  and  natural  space  boundaries  are  determined  from  latitude- 
longitude  coordinates  contained  in  an  external  file  supplied  with  the  program. 
After  the  user  enters  the  northwest  corner,  the  geographic  location  of  the  grid 
is  determined  and  drawn.  Station  location  symbols  are  then  plotted  on  the  map. 
This  map  forms  the  background  for  all  the  contoured  fields  generated  by  the  model. 

The  northwest  corner  grid  coordinates  are  then  passed  to  linked  programs  that 
calculate  the  station’s  location  relative  to  the  grid  and  then  determines  the  radii  of 
the  ten  closest  stations  from  each  grid  point.  Stations  are  selected  automatically 
from  a  list  of  North  American  rawinsonde  stations,  and  those  selected  are  located 
inside  of  or  within  one  mesh  length  outside  the  boundary  of  the  18  by  15  grid. 

The  ten  smallest  radii  to  stations  at  each  grid  point  are  then  passed  to  another 
program  which  determines  the  influence  of  each  station’s  observation  according  to 


Since  the  format  of  the  required  data  varies  depending  upon  the  source  (i.e.,  WMO 
code,  files  from  the  NMC,  or  some  other  source),  the  Barnes’  model  provides  three 
different  data  ingest  routines. 

The  first  method  is  to  manually  enter  all  the  data.  The  second  mode  allows 
input  from  the  NMC  mandatory  level  plot  files  provided  that  they  can  be  made 
available  on  a  5  1  /4  inch  floppy  diskette  for  the  IBM  series  of  microcomputers. 
The  final  mode  uses  files  that  reside  on  computers  which  monitor  the  FAA  604  line 
and/or  acquire  sounding  data  through  a  link  with  the  PROFS  data  ingest  facility. 
The  data  in  the  third  mode  are  available  through  the  use  of  a  1200  baud  modem 
over  a  telephone  link. 

A  considerably  more  difficult  problem  is  data  coverage  and  completeness. 
The  grid  point  interpolation  scheme  cannot  be  relied  upon  to  give  meaningful 
information  either  at  the  boundaries  of  the  data  domain  or  in  areas  missing  data, 
which  represent  internal  boundaries.  The  scheme  is  most  reliable  when  it  has  a 
quasi-uniform  distribution  of  data  over  the  grid  ;  d  data  outside  the  boundaries 
of  the  grid.  Barnes  (1986b)  pointed  out  that  a  minimum  of  about  75%  of  the 
observations  over  the  data  domain  are  required  in  order  to  obtain  reliable  diagnostic 
results. 

One  method  for  enhancing  the  data  coverage  over  oceans  adjacent  to  North 
America  is  to  incorporate  “bogus”  data  points  at  selected  latitude/longitude  points. 
Representative  data  for  the  bogus  points  can  be  obtained  from  contoured  mandatory 
pressure  level  maps  as  provided  by  the  NMC  on  AFOS.  The  accuracy  of  the  bogus 
points  needs  to  be  no  greater  than  —  10  m.  The  data  ingest  procedures  would  have 
to  be  modified  slightly  to  incorporate  these  bogus  points.  This  was  not  done  in  the 
original  program  since  the  number  of  bogus  points  used  would  be  dependent  on  the 
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selected  geographical  area  of  the  grid. 

Another  solution  to  the  data  coverage  and  accuracy  problem  is  to  use  the  NMC 
grid  point  values  obtained  from  LFM  initial  conditions  or  predicted  height  fields. 
For  more  information  on  this  topic  see  Barnes  (1986a). 

For  manual  data  entry  the  user  is  prompted  (by  messages  on  the  screen) 
for  the  required  data  station  by  station.  The  order  in  which  data  are  entered  is 
important  since  the  interpolation  scheme  assumes  that  the  first  element  in  a  data 
file  is  associated  with  station  “x”,  the  second  element  is  for  station  “y”,  and  so  on. 

Format  of  the  data  input  is  the  same  as  it  appears  on  the  NMC  mandatory 
pressure  charts.  Three  digits  are  entered  for  the  height  while  temperature  and  dew 
point  depression  are  entered  as  they  appear  on  the  charts.  If  the  WMO  coded  data 
are  used,  then  the  tenths  of  values  for  temperature  and  dew  point  depression  should 
also  be  entered  for  greater  accuracy  in  subsequent  lapse  rate  computations. 

Missing  data  are  entered  as  four  nines  (9999)  for  each  observation  that  is  not 
available  or  is  judged  to  be  in  error.  After  all  the  data  have  been  entered,  the  editing 
phase  of  the  program  is  entered.  Subprogram  “EDITUAIR”  allows  corrections  of 
any  of  the  data.  In  the  editing  program  instructions  are  given  as  to  how  to  flag  and 
correct  errors.  All  the  data  entered  are  redisplayed  on  the  screen  for  the  user  to 
verify  and  correct  if  necessary.  Data  must  be  entered  as  four  digits  for  the  height 
fields  in  the  editing  phase  and  not  the  three  digit  number  used  in  the  manual  entry 
portion  of  the  program. 

As  the  data  are  displayed  on  the  screen,  an  arrow7  may  appear  next  to  a  value. 
This  signifies  that  the  item  is  outside  predetermined  limits.  These  limits  are  stored 
in  an  external  file  and  the  file  can  be  edited  by  the  user  to  change  these  limits  at 
any  time  before  running  the  program. 


Upon  completion  of  the  editing  (correcting)  phase,  the  program  is  automatic 
and  there  isn’t  any  more  user  interaction  until  a  menu  for  map  display  is  presented. 
The  program  is  run  by  calling  the  required  subprograms. 

The  subprogram  “LapseRate”  computes  the  lapse  rate  of  temperature  with 
pressure  between  the  four  mandatory  levels  based  upon  the  input  temperature  data 
and  the  computed  thickness  of  each  layer.  If  the  values  are  entered  to  the  nearest 
degree,  as  are  usually  obtained  from  charts,  the  lapse  rates  can  exceed  dry  adiabatic 
(9.8°C/km)  on  some  occasions.  This  is  due  to  the  truncation  of  tenths  values  in  the 
input  data.  This  will  cause  the  magnitudes  of  lapse  rates  as  seen  on  the  contoured 
charts  to  be  in  error  somewhat,  but  the  relative  distribution  of  lapse  rate  values  can 
be  believed.  One  can  still  see  where  there  are  centers  of  maximum  and  minimum 
values. 

Since  there  are  only  three  layers  in  which  lapse  rates  can  be  calculated,  and 
because  the  interpolation  programs,  for  symmetry  reasons,  are  designed  to  handle 
four  variables  at  a  time,  this  program  will  file  the  850  mb  temperature  in  addition 
to  the  three  lapse  rate  fields. 

Subprogram  “MixRatio”  computes  the  mixing  ratio  fields  for  the  four  pressure 
levels  previously  defined  while  subprogram  “InterpZ”  interpolates  geopotential 
heights  to  grid  points.  The  interpolation  scheme  used  is  a  two-pass  Barnes  analysis 
for  which  the  weight  function  parameter  (C  =  63000  km2  as  in  Barnes,  1985)  is  not 
changed  for  the  second  pass.  This  is  a  departure  from  the  weighting  scheme  used 
by  Barnes  (1985)  which  was  determined  to  analyze  the  scales  of  height  variations 
that  were  very  close  to  the  minimum  detectable  signals  obtainable  from  the  U.S. 
rawinsonde  network. 

The  program  also  interpolates  the  temperature  lapse  rate  (“InterpLR”)  and 
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mixing-ratios  (“InterpMR”).  Then  “ThicknMZ”  computes  the  thickness  and  mean 
height  for  each  of  the  three  layers. 
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Subprogram  “QVector”  computes  for  each  of  the  three  layers  the  thickness 
temperature,  geostrophic  temperature  advection,  Q-vector  components,  divergence 
of  Q,  and  the  geostrophic  frontogenesis  function.  Each  of  these  diagnostic  fields  is 
computed  from  various  combinations  of  gradients  of  geostrophic  wind  at  the  mean 
height  of  the  layer  and  layer  temperature.  Layer  temperature  values  are  computed 
from  thickness  values  for  each  layer.  All  gradient  values  are  subjected  to  a  nine-point 
smoothing  routine  before  entering  subsequent  computations. 

Subprogram  “Stability”  uses  the  mean  potential  temperature  of  each  layer 
to  compute  the  stability  of  each  layer  expressed  as  a  rate  of  change  of  potential 
temperature  with  decreasing  pressure.  Thus,  the  fields  should  always  be  positive. 
Negative  values  would  indicate  greater  than  adiabatic  lapse  rates. 

Subprogram  “Vorticity”  computes  for  each  of  the  three  layers  the  absolute 
geostrophic  vorticity.  The  computed  geopotential  Laplacians  are  subjected  to 
a  nine-point  smoothing  routine.  The  vorticity  fields  are  computed  from  the 
mean  height  fields  for  the  three  layers  which  have  mean  pressure  values  of  771.4, 
591.6,  and  387.3  mb.  The  mean  pressure  values  are  obtained  with  the  equation 
MP  =  EX P((log(P2)+log(Pl))/2).  In  this  equation  MP  stands  for  mean  pressure, 
PI  is  the  lower  (in  height)  pressure  value,  and  P2  is  the  upper  pressure  level.  Using 
an  arithmetic  mean  would  have  been  more  appropriate  since  pressure  coordinates 
were  used  in  the  program.  However,  in  this  case,  the  mean  pressure  values  are  only 
slightly  different  from  what  would  be  obtained  using  the  arithmetic  mean  (775,  600, 
and  400  mb).  For  convenience,  these  and  other  fields  based  upon  mean  height  values 
carry  labels  of  800,  600,  and  400  mb  respectively. 
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Both  “SWLower”  and  “SWUpper”  estimate  the  local  lapse  rate  (stability) 
tendency  in  a  layer.  The  lower  layer  is  from  850  to  500  mb  and  the  upper 
layer  is  from  700  to  300  mb.  This  overlap  arises  because  the  omega  forcing  is 
obtained  at  three  levels.  The  forcing  at  600  mb  is  used  in  both  the  lower  and 
upper  tendency  computations.  The  computation  is  based  upon  the  expression  for 
adiabatic  processes  and  only  includes  the  effect  of  horizontal  advection  of  stability 
by  the  geostrop  hie  wind. 

The  next  subprogram  “StabTend”  computes  the  stability  tendency  and  files  it 
in  two  forms.  One  is  in  the  dimensions  in  which  it  is  computed  and  the  other  is  in 
dimensions  more  familiar  to  operational  meteorologists,  that  is,  °C  km-1  (12h)"“1. 

The  lapse  rate  fields  themselves  are  based  upon  reported  temperature  data  at 
mandatory  levels.  The  computed  lapse  rate  tendency  is  based  upon  height  data 
alone,  through  application  of  quasi-geostrophic  concepts.  Thus,  the  changes  in 
lapse  rate  fields  from  one  time  to  the  next  will  not  necessarily  be  similar  to  the 
computed  tendency  fields.  It  is  unknowm  at  this  time  whether  their  differences  will 
reveal  anything  about  the  diabatic  and  ageostrophic  advective  processes  that  are 
not  included  in  quasi-geostrophic  theory. 

The  next  two  programs  “GTLower”  and  “GTUpper”  input  the  fields  required 
to  compute  geopotential  tendencies  in  the  subprogram  “GeopTend.”  The  tendencies 
are  computed  as  a  scaled  advection  of  vorticity. 

This  program  also  computes  two  forecast  fields.  These  are  height  fields  for  t+12 
hours  and  are  computed  from  the  initial  height  fields  at  700  and  500  mb  and  the 
tendency  fields  for  the  upper  and  lower  layers.  The  forecast  height  fields  represent 
a  gross  estimate  of  how  the  height  fields  would  change  over  the  next  few  hours  if 
geostrophic  forcing  were  the  only  influence.  The  forecast  height  fields  are  the  sum 


of  the  initial  height  fields  and  the  geopotential  tendencies  discussed  above. 

The  remainder  of  the  programs  contained  within  the  batch  file  deals  with 
the  graphics  output.  The  package  provided  by  Dr.  Barnes  contains  two  display 
programs.  As  previously  mentioned  twenty  maps  are  automatically  contoured  on 
the  background  map.  Once  the  maps  have  been  generated  they  can  be  individually 
retrieved  through  the  menu  described  below,  and  they  can  be  output  to  a  printer. 
Additionally,  twenty  other  maps  can  be  contoured  from  an  alternate  menu.  These 
additional  maps  can  also  be  sent  to  the  printer. 

Program  “ContPkg”  (Contour  Package)  contains  the  filenames  of  maps  to  be 
contoured  in  DATA  statements.  The  contour  intervals  of  each  map  are  contained 
in  a  file  named  “Contour.DAT”  and  can  be  modified  by  the  user  before  running  the 
program. 

The  present  program  will  blow  up  if  too  many  contours  are  present.  This  is 
due  to  the  limited  array  space  available  for  contour  points.  In  order  to  alleviate  this 
problem  another  program  has  been  added  that  checks  the  maximum  and  minimum 
values  to  be  contoured  and  then  ensures  that  the  program  won’t  try  to  print  too 
many  contours.  If  the  default  contour  values  are  such  that  more  than  ten  different 
contours  values  will  result,  then  the  contour  interval  is  adjusted  automatically.  The 
number  of  different  contour  values  to  be  drawn  can  be  changed  by  adjusting  one 
parameter  in  the  program  “AdjCntr”  (Adjust  Contour). 

“Map View”  provides  the  main  map  viewing  menu  for  the  20  preselected  fields. 
The  individual  maps  can  be  displayed  in  any  order  and  also  “tagged”  for  later  output 
to  the  printer.  There  is  also  an  option  that  will  automatically  print  all  20  maps. 
Since  the  maps  have  already  been  contoured  and  saved  to  disk  the  regeneration  of 
these  maps  is  almost  instantaneous. 


“AltMenu”  is  another  map  viewing  program  that  can  be  selected  from  the  main 
menu  described  above.  The  difference  with  the  above  program  is  that  these  20  maps 
are  contoured  while  the  user  watches  instead  of  being  retrieved  from  disk. 

Since  other  maps  may  be  desired  by  the  user,  AltMenu  can  be  modified  to  give 
the  user  more  maps  from  which  to  choose.  This  could  be  done  by  chaining  another 
menu  to  the  program.  The  only  difference  in  this  new  menu  program  would  be 
the  contoured  maps  that  the  user  could  select  from.  The  DATA  statements  that 
contain  the  data  files  to  be  contoured  would  have  to  be  changed.  This  could  be 
done  as  many  times  as  desired,  so  that  all  the  data  generated  could  be  contoured. 

The  programs  described  above  used  centered  difference  forms  to  compute 
derivatives  at  all  interior  points.  Boundary  points  were  evaluated  by  extrapolating 
outward  from  the  two  adjacent  interior  points.  The  four  corner  points  used  the 
two  values  on  a  diagonal  line  with  the  respective  corner.  The  remainder  of  the 
boundary  points  used  values  either  on  a  horizontal  or  vertical  line  depending  on  the 
point  being  evaluated. 

The  Barnes’  program  provides  a  picture  of  vertical  motion  at  three  levels  in 
the  atmosphere  while  the  standard  NGM  facsimiles  provide  vertical  motion  only 
for  the  700  mb  level.  Having  a  picture  of  vertical  motion  at  three  levels  gives  the 
meteorologist  a  more  comprehensive  picture  of  the  state  of  the  atmosphere.  For 
instance,  if  there  was  low  level  upward  motion  and  convection  was  possible,  upper 
level  subsidence  may  prevent  or  may  cap  thunderstorm  development. 

The  Barnes’  program  is  easily  run  and  it  is  fast.  The  program  can  be  run  on  a 
microcomputer  which  is  usually  available  in  any  weather  station  today.  After  data 
ingest,  the  user  can  have  diagnostic  charts  within  10  min.  Since  most  upper  air 
data  are  available  between  2  and  3  h  after  data  time,  the  user  can  have  the  charts 
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available  at  least  1  h  before  the  commonly  used  charts  are  available.  It  is  usually  4 
h  after  data  time  that  the  first  of  the  diagnostic  or  prognostic  products  are  available 
in  chart  form. 

Very  little  user  interaction  is  required  after  data  ingest.  This  makes  it  more 
likely  that  the  operational  meteorologist  will  use  the  program.  Most  of  the  user 
interaction  required  is  obtaining  and  verification  of  the  input  data.  The  Barnes’ 
program  has  defined  limits  on  the  values  of  input  data  and  flags  any  data  outside 
of  these  limits.  This  makes  the  data  verification  a  much  easier  process. 

All  computations  in  the  Barnes  program  used  pressure  as  the  vertical  coordi¬ 
nate.  The  derivation  section  presented  earlier  is  using  a  modified  x,  y,  z  coordinate 
system.  In  the  x,  y,  p  coordinate  system  (12)  becomes 

V2(<ruO  +  /2|^  =  -2(V.Q) 

where  <r  =  —  {a  /  Q)(d6  /  dp)  and 


dx  \dp  )  ’  dy  V  dp 
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CHAPTER  IV 

CASE  STUDIES 


THUNDERSTORMS  IN  TEXAS  ON  17  MARCH  1987 


On  17  March  1987,  eastern  Texas  experienced  several  tornadoes,  hail,  and 
high  winds.  The  most  severe  weather  produced  by  the  storm  system  occurred 
between  1100  UTC  and  1400  UTC  (0500  and  0800  LST).  This  case  study  will  look 
at  a  comparison  of  how  the  NWS  prognostic  maps  compare  with  the  results  of  the 
Barnes’  3-Level  Q- vector  program.  First,  a  discussion  of  the  synoptic  situation  will 
be  presented. 

The  1200  UTC  surface  map  (not  shown)  showed  a  deepening  992  mb  low 
approximately  80  km  northeast  of  Amarillo,  Texas  with  a  warm  front  extending 
eastward  through  central  Oklahoma  and  northern  Arkansas.  A  cold  front  extended 
from  the  low  southeastward  to  the  Austin,  Texas  area  and  then  southward  to 
Brownsville,  Texas.  The  surface  map  analysis  also  showed  a  pressure  trough  or 
possible  convergence  line,  to  the  east  of  the  cold  front,  that  extended  southward 
from  the  Wichita  Falls,  Texas  area  to  Tyler,  Texas  and  southwestward  from  Tyler 
to  about  40  km  west  of  Houston.  The  severe  weather  occurred  in  the  vicinity  of 
this  trough  line. 

The  NWS  Radar  Summary  at  1035  UTC  (not  shown)  showed  a  severe  weather 
watch  area  that  coincided  with  the  area  of  the  pressure  trough.  The  severe  weather 
watch  had  been  issued  6  h  earlier  and  was  due  to  expire  at  1100  UTC.  The  summary 
also  showed  thunderstorms  throughout  eastern  Texas  with  maximum  cloud  tops 
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reaching  to  45,000  ft.  Subsequent  radar  summaries  showed  these  thunderstorms 
moving  eastward  and  the  heights  of  maximum  cloud  tops  increasing. 

The  upper  air  maps  showed  strengthening  of  the  low  pressure  system  in  the 
lower  and  upper  layers  while  weakening  in  the  middle  layers.  The  850  mb  1200 
UTC  chart  showed  a  minimum  height  of  1332  gpm  near  Amarillo,  Texas  which  was 
30  m  lower  than  on  the  0000  UTC  chart  on  the  same  day.  There  was  strong  cold 
advection  in  the  rear  of  the  low.  The  700  mb  chart  showed  weak  cold  advection  in 
the  southern  part  of  the  low  while  the  500  mb  chart  showed  mostly  warm  advection 
in  the  rear  of  the  low  along  with  an  area  of  colder  air  to  the  southwest  side  of  low 
center.  The  300  mb  chart  had  cold  advection  into  the  northwestern  side  of  the  low. 
Figures  6  through  9  show  these  charts. 


FIG.  6.  850  mb  height  and  temperature  for  1200  UTC  17  March  1987.  Height 
(solid  lines)  is  in  dam  while  temperature  (dashed  lines)  is  in  °C. 


Concerning  these  and  subsequent  figures,  contour  values  are  only  labeled  on 
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FIG.  7.  700  mb  height  and  temperature  for  1200  UTC  17  March  198/ 

Labelling  is  the  same  as  in  Fig.  6. 
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FIG.  9.  300  mb  height  and  temperature  for  1200  UTC  17  March  1987. 

Labelling  is  the  same  as  in  Fig.  6. 


the  edges  of  the  charts.  Maximum  and  minimum  locations  have  different  markings 
depending  upon  the  sign  of  each  value.  If  the  two  values  are  of  opposite  signs 
then  maxima  are  labeled  with  an  and  minima  are  labeled  with  a  while 
values  with  the  same  sign  result  in  a  “X”  and  “N”  respectively.  Any  maximum  or 
minimum  values  mentioned  in  the  following  discussion  are  obtained  from  separate 
numerical  data  files  that  list  these  values  along  with  their  grid  point  locations. 

The  NGM  36  hour  forecast  vorticity  chart  (valid  1200  UTC  17  March  1987,  not 
shown),  based  upon  the  0000  UTC  16  March  data  showed  a  vorticity  maximum  in 
central  Texas  with  weak  positive  vorticity  advection  (PVA)  extending  into  northern 
central  Texas  and  negative  vorticity  advection  (NVA)  into  eastern  Texas.  On  the 
subsequent  NGM  (12  h  later),  the  24  h  forecast  vorticity  chart  showed  the  maximum 
in  about  the  same  location  but  the  magnitude  of  the  vorticity  values  was  higher.  The 
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vorticity  advection  in  eastern  Texas  was  near  zero.  The  NGM  chart,  computed  from 
the  0000  UTC  17  March  data,  was  not  available  due  to  communication  problems. 
This  would  have  produced  a  forecast  vorticity  chart  computed  from  the  latest 
available  data  before  the  onset  of  severe  weather. 

The  NGM  vorticity  charts  are  combined  with  the  500  mb  height  field  forecast 
for  the  same  period.  The  Barnes’  3- Level  Q-vector  program  (hereafter  referred  to  as 
the  Barnes  program)  also  computes  vorticity  but  the  contoured  fields  are  combined 
with  thickness  values  for  the  layer.  The  Barnes  program  computes  vorticity  for  the 
600  mb  level  and  computes  the  thickness  based  upon  the  700  mb  and  500  mb  data. 
Trenberth  (1978)  presented  arguments  for  evaluating  vertical  motion  based  upon 
advection  of  vorticity  by  the  thermal  wind.  Since  the  thermal  wind  is  proportional 
to  thickness,  the  user  may  choose  to  evaluate  the  vertical  motion  field  using  this 
vorticity  chart.  In  addition  to  the  600  mb  level,  vorticity  charts  are  also  generated 
for  the  800  mb  and  400  mb  levels  with  the  appropriate  thickness  contours. 

Figures  10  and  11  show  the  600  mb  vorticity  field  combined  with  the  thickness 
field.  Figure  10  is  for  0000  UTC  17  March  data  and  shows  that  there  is  PVA  across 
northeastern  Texas  and  into  Lousiania.  Figure  11  is  for  the  1200  UTC  data  and 
shows  that  the  PVA  across  northeastern  Texas  has  become  very  weak.  The  Barnes 
vorticity  charts  differ  from  the  NGM  charts  as  to  a  vorticity  maxima  location  and 
in  advection  thro  gh  Texas.  Part  of  the  discrepancy  may  be  due  to  the  NGM  wind 
observations  and  predictions  being  subjected  to  smoothing  algorithms. 

Figures  12  and  13  show  the  diagnostic  geopotential  tendency  fields  at  the  500 
mb  and  the  700  mb  level,  measured  in  decameters  per  12  h,  at  0000  UTC.  Figures 
14  and  15  show  the  same  fields  12  h  later.  Comparing  these  charts  shows  that,  at 
the  700  mb  level,  the  magnitude  of  the  falls  decreased  over  most  of  Texas  over  the 
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FIG.  10.  600  mb  vorticity  and  700-500  thickness  for  0000  UTC  17  March 

1987.  Vorticity  isopleths  are  dashed  and  the  thickness  (dam)  are  solid  lines. 


FIG.  11.  600  mb  vorticity  and  700-500  thickness  for  1200  UTC  17  March 

1987.  Labelling  is  the  same  as  in  Fig.  10. 
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period  from  0000  UTC  to  1200  UTC.  The  500  mb  charts  show  eastern  Texas  under 
j  falling  tendency  at  0000  UTC  and  rising  tendency  12  h  later  except  in  extreme 

eastern  Texas.  The  500  mb  geopotential  tendency  through  Texas  and  Lousiania  is 
in  aggreement  with  the  positive  vorticity  advection  (PVA)  pattern  shown  in  Fig. 
|  11.  The  areas  of  PVA  coincide  with  the  areas  of  falling  geopotential  tendency.  The 

700  mb  negative  tendency  shows  that  the  lower  layers  were  having  geopotential  falls 
through  Texas  while  the  500  mb  charts  show  that  there  were  geopotential  rises  in 
j  the  middle  layers  over  western  Texas.  The  Barnes’  model  also  computes  stability 

|  tendency.  The  charts  for  the  700  mb  and  500  mb  layers  (not  showm)  show  that  the 

i 

i  700  mb  level  was  becoming  more  unstable  and  the  500  mb  level  was  becoming  more 

i  stable. 


FIG.  12.  500  mb  geopotential  tendency  for  0000  UTC  17  March  1987.  Units 

are  dam  (12  h)_1 

The  V  •  Q  charts  for  the  800  mb,  600  mb,  and  400  mb  layers  (Figs.  16-18) 


FIG.  15.  700  mb  geopotential  tendency  for  1200  UTC  17  March  1987. 

Labelling  is  the  same  as  in  Fig.  12. 


show  the  vertical  motion  associated  with  the  storm  system  at  1200  UTC.  The  zero 
contour  line  on  the  600  mb  chart  (Fig.  17)  fairly  well  coincides  with  the  surface  cold 
front  position.  All  three  charts  show  relatively  strong  subsidence  (positive  values) 
over  the  western  two  thirds  of  Texas  and  rising  motion  (negative  values)  over  the 
remainder  of  Texas.  The  rising  motion  over  the  eastern  part  of  Texas,  in  the  mid 
troposphere,  is  also  consistent  with  the  idea  presented  by  Djuric  (1969)  that  upward 
vertical  motion  coincides  with  PVA  above  the  level  and  strong  warm  air  advection 
below  the  level.  The  rising  motion  through  Oklahoma,  Kansas,  and  Nebraska  can 
be  associated  with  an  inverted  trough,  extending  through  the  center  of  these  states, 
at  the  surface. 

Water  vapor,  another  important  variable  in  a  storm  environment,  is  depicted 
as  mixing  ratio  in  Figs.  19-24.  Comparing  the  850  mb  mixing  ratio  at  0000  UTC 
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FIG.  18.  400  mb  V  •  Q  for  1200  UTC  17  March  1987.  Labelling  is  the  same 

as  in  Fig.  16. 


(Fig.  19)  with  1200  UTC  (Fig.  20)  shows  a  increasingly  humid  and  pronounced 
ridge  over  eastern  Texas  and  Lousiania  extending  up  through  the  central  plains 
states.  The  water  vapor  patterns  show  the  same  results  when  comparing  the  700 
mb  charts  (Figs.  21  and  22)  and  the  500  mb  charts  (Figs.  23  and  24).  In  western 
Texas,  on  all  the  charts,  the  amount  of  moisture  had  decreased  over  the  12  h  period. 

The  subsidence,  shown  on  the  V  •  Q  charts,  and  the  decreasing  humidity 
over  western  Texas  showed  as  rapid  clearing  behind  the  front.  The  36  h  NGM 
precipitation  and  vertical  velocity  forecast  charts,  valid  at  1200  UTC  17  March 
1987  showed  just  the  opposite.  They  showed  continuing  precipitation  and  lifting 
motion.  The  48  h  charts  computed  from  the  same  data  base  showed  western  Texas 
under  subsidence  and  decreased  precipitation.  The  subsequent  precipitation  chart 
shows  the  precipitation  still  occurring  at  1200  UTC  and  ending  12  h  later.  Since  the 
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FIG.  19.  850  mb  mixing  ratio  for  0000  UTC  17  March  1987.  Units  are  g 
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FIG.  20.  850  mb  mixing  ratio  for  1200  UTC  17  March  1987.  Labelling  is  the 
same  as  in  Fig.  19. 
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FIG.  21.  700  mb  mixing  ratio  for  0000  UTC  17  March  1987.  Labelling  is  the 
same  as  in  Fig.  19. 
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FIG.  22.  700  mb  mixing  ratio  for  1200  UTC  17  March  1987.  Labelling  is  the 
same  as  in  Fig.  19. 
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NGM  charts  computed  from  the  0000  UTC  17  March  1987  data  was  unavailable, 
it  it>  unknown  if  the  NGM  prognostic  maps  were  still  approximately  12  h  slow. 
Based  upon  continuity  of  the  NGM  products  it  is  believed  that  the  NGM  products 
continued  to  be  slow  in  their  timing. 

The  forecaster  using  the  NGM  products  computed  from  the  0000  UTC  17 
March  data,  for  example,  could  have  compared  these  products  to  the  Barnes  model 
results.  If,  as  assumed,  the  NGM  continued  the  precipitation  over  western  Texas, 
then  the  forecaster  could  have  adjusted  his  use  of  the  NGM  products  based  upon 
the  quasi-geostrophic  forcing  as  depicted  by  the  Barnes  program  and  comparison 
with  available  satellite  pictures.  The  forecaster  would  have  seen  the  decreasing 
precipitation  over  western  Texas  and  the  enhanced  storm  environment  over  eastern 
Texas.  The  rising  motion  and  increased  humidity  ahead  of  the  front  helped  enhance 
the  storm  environment  in  that  area. 

From  the  Q  and  temperature  charts  (Figs.  25-27)  the  V  ■  Q  maps,  previously 
discussed,  can  be  constructed.  The  Barnes’  program  labels  the  charts  with 
part  of  the  title  as  ’’Thickness  Temperature.”  The  contoured  field  is  a  computed 
temperature  based  upon  the  thickness  of  the  layer.  There  is  rising  motion  where 
the  Q-vectors  converge  and  conversly,  subsidence  where  they  diverge. 

Figure  26  shows  the  convergence  of  the  Q-vectors  over  eastern  Texas  and 
western  Lousiania.  This  corresponds  to  the  rising  motion  depicted  in  the  same  area 
in  Fig.  17.  Western  Texas  shows  V  •  Q  >  0  which  corresponds  to  the  subsidence 
area  in  Fig.  17. 

Comparing  the  direction  of  the  Q-vector  with  respect  to  the  temperature 
contours  yields  the  GFF.  Wfien  the  Q-vector  has  a  component  pointing  from  colder 
air  toward  warmer  air,  the  potential  temperature  gradient  is  enhanced  and  a  positive 
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FIG.  25.  800  mb  Q  and  temperature  for  1200  UTC  17  March  1987.  The 

temperature,  computed  from  thickness  of  the  layer  from  850  mb  to  700  mb,  is  in 
°C.  Arrows  represent  the  direction  and  magnitude  of  Q. 

GFF  results.  A  positive  GFF  signifies  frontogenesis.  A  component  of  Q  from  warm 
to  cold  air  indicates  a  geostrophic  frontolytic  situation  (negative  GFF). 

Figures  28-30  show  the  resulting  GFF  for  the  800  mb,  600  mb,  and  400  mb 
levels.  Comparing  Fig.  25  with  Fig.  28  shows  that  across  central  Texas  at  800  mb 
the  Q-vectors  point  from  cold  to  warm  and  have  a  corresponding  positive  GFF.  In 
this  case  the  motion  was  tending  to  intensify  the  thermal  gradient  in  the  vicinity 
of  the  frontal  system.  A  comparison  of  Fig.  26  with  Fig.  29  shows  that  at  600  mb 
there  is  a  weak  component  of  motion  from  cold  to  warm  across  central  Texas  and 
therefore  a  relatively  weak  positive  area  of  GFF  in  the  same  region.  At  the  400  mb 
level  the  northeastern  part  of  Texas  has  a  positive  GFF. 

One  must  be  careful  when  interpreting  the  GFF.  In  this  situation  the  GFF  is 
acting  to  maintain  the  frontal  circulations  and  not  to  support  the  actual  formation  of 
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FIG.  26.  600  mb  Q  and  temperature  for  1200  UTC  17  March  1987.  Labelling 
is  the  same  as  in  Fig.  25. 
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FIG.  27.  400  mb  Q  and  temperature  for  1200  UTC  17  March  1987.  Labelling 
is  the  same  as  in  Fig.  25 
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FIG.  30.  400  mb  GFF  for  1200  UTC  17  March  1987.  Labelling  is  the  same 

as  in  Fig.  28. 

a  front.  As  Holton  (1979  pp.  240-242)  pointed  out,  the  quasi-geostrophic  equations 
are  invalid  beyond  the  initial  stages  of  frontogenesis.  Barnes  (1985)  hypothesized 
that  geostrophic  frontogenesis  might  provide  the  organizing  environment  in  which 
the  ageostrophic  frontogenetic  components  operate  more  effectively. 

The  charts  in  this  case  show  that  the  front  was  intensifying  in  the  lower  layers 
and  that  there  was  rising  motion  at  each  level  over  eastern  Texas.  The  strong 
vertical  motion  combined  with  the  influx  of  water  vapor  allowed  the  thunderstorms 
to  produce  the  severe  weather. 

Now  that  it  has  been  shown  how  the  Barnes’  diagnostic  charts  corresponded 
with  the  prognostic  charts,  one  may  look  closer  at  how  the  diagnostic  charts  related 
to  the  actual  synoptic  conditions  at  the  time.  This  will  consist  of  comparing  the 
V  •  Q  charts  with  the  NWS  Radar  Summaries.  To  start  this  procedure,  a  chart 


can  be  used  with  the  appropriate  segment  of  the  U.S.  outlined  and  marked  off  in 
1°  latitude/longitude  increments.  The  NWS  R.adar  Summary  can  be  overlaid  onto 
this  chart.  The  only  concern  is  with  the  outline  of  the  precipitation  area.  The 
last  preliminary  step  is  to  overlay  the  V  •  Q  charts  for  the  800,  600,  and  400  mb 
levels  onto  the  same  chart.  Unless  otherwise  specified  in  the  following  discussion, 
the  contour  drawn  is  the  V  •  Q  =  0  contour.  The  overlaid  chart,  constructed  for 
1200  UTC  17  March  1987,  is  shown  in  Fig.  31.  For  the  sake  of  clarity,  the  1° 
latitude/longitude  increments  are  omitted  in  Fig.  31. 
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FIG.  31.  Composite  chart  for  1200  UTC  17  March  1987.  Precipitation  (shaded 
area)  is  overlaid  with  V  •  Q  computed  at  800  mb  (short  dashed  lines),  600  mb  (long 
dashed  lines),  and  400  mb  (dotted  lines).  The  arrows  represent  the  side  of  the 
contour  with  V  •  Q  <  0. 

A  comparison  will  be  made  of  implied  vertical  motion  with  the  precipitation 
areas  by  counting  the  latitude/longitude  points  that  are  within  regions  meeting  the 
desired  criteria.  Two  sets  of  comparisons  will  be  made  using  the  composite  map. 
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The  first  set  of  comparisons  determines  the  ratio  of  the  number  of  points  with 
precipitation  and  lifting  to  the  total  number  of  precipitation  points.  This  is  done 
for  the  800,  600.  and  400  mb  levels.  Also,  a  ratio  of  the  number  of  points  with 
precipitation  and  lifting  at  any  level  to  the  total  number  of  precipitation  points  will 
be  determined.  For  convience  this  final  ratio  will  be  called  a  composite  ratio. 

The  second  set  of  comparisons  determines  the  ratio  of  the  number  of  implied 
subsidence  points  that  are  precipitation  free  to  the  total  number  of  implied  subsi¬ 
dence  points.  As  above,  this  will  be  done  for  each  of  the  three  levels.  Finally  a 
composite  ratio  will  be  determined  by  comparing  the  precipitation  free  subsidence 
points  to  number  of  points  with  implied  subsidence  at  every  level. 

Table  1  shows  the  results  of  the  two  sets  of  comparisons.  The  first  line  is  for 
the  comparisons  dealing  with  precipitation  and  lifting.  The  second  fine  is  for  the 
precipitation  free  areas  and  subsidence  comparisons.  There  are  two  numbers  in  each 
column  of  both  lines.  The  first  number  represents  the  number  of  points  meeting 
the  required  criteria  and  the  second  number  represents  the  total  number  of  points 
considered. 

Looking  at  Table  1  in  more  detail,  it  is  seen  in  the  precipitation  case  that  at 
800  mb,  68  out  of  85  precipitation  points  are  in  a  area  of  implied  lifting.  Rounding 
this  to  the  nearest  percentage,  this  is  80%.  The  percentages  for  each  of  the  ratios 
in  Table  1  is  given  in  Table  2. 

It  is  seen,  from  Tables  1  and  2,  that  the  results  for  each  level  were  almost 
the  same.  It  is  also  seen,  from  the  composite  map  (Fig.  31),  that  this  is  a  fairly 
homogeneous  area  of  precipitation.  Why  weren’t  the  figures  in  Tables  1  and  2 
higher?  Part  of  the  reason  may  be  due  to  the  smoothing  and  interpolation  schemes 
used  in  the  evaluation  of  the  V  •  Q  charts. 


TABLE  1.  Ratios  of  precipitation  points  with  lifting  over  total  precipitation 
points  and  the  ratios  of  precipitation  free  subsidence  points  over  total  subsidence 
points  for  1200  UTC  17  March  1987.  For  each  group  of  two  numbers  the  first  number 
represents  the  number  of  points  meeting  the  required  criteria  and  the  second  number 
is  the  number  of  to^al  points  considered.  The  composite  column  considers  points 
with  lifting  at  any  level. 
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TABLE  2.  Percentages  for  data  in  TABLE  1. 


800  mb  600  mb  400  mb  Composite 


78  91 

91  95 


600  mb 


400  mb 


An  additional  test  was  performed  on  the  precipitation  case  for  800  mb  in  order 
to  determine  if  the  percentage  shown  in  Table  2  was  statistically  significant.  The 
confidence  interval  (99%)  for  a  binomial  experiment  was  found  to  be  <  60%.  The 
result  shown  in  Table  2  is  that  80%  of  the  precipitation  areas  had  upward  motion 
at  the  800  mb  level  and  therefore  the  result  is  statistically  significant.  This  was 
expected  since  low  level  upward  motion  is  expected  in  precipitation  regions. 

The  second  line,  in  Tables  1  and  2,  shows  a  higher  percentage  of  subsidence 
areas  being  precipitation  free.  Since,  in  this  case,  the  precipitation  points  within 
subsidence  areas  were  near  the  V  •  Q  =  0  contour,  one  may  possibly  infer  that  the 
precipitation  was  abating  in  these  areas.  For  stronger  magnitudes  of  subsidence 
(i.e.  V  •  Q  >  4  *  10-3 '  s-3mi-1 )  the  subsidence  areas  are  precipitation  free.  To  see 
this  compare  the  values  of  V  •  Q  in  Figs.  16-18  with  the  precipitation  area  in  Fig. 
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Another  set  of  comparisons  was  made  in  order  to  ascertain  the  relationship 
between  values  of  V  •  Q  and  precipitation.  The  first  comparison  was  to  see  what 
proportion  of  the  upward  vertical  motion  area  had  precipitation.  Once  again,  the 
V  •  Q  =  0  contour  is  used. 

The  second  comparison  was  as  above,  except  with  use  of  a  stronger  value  of 
vertical  motion,  i.e.  a  critical  value  of  V  •  Q  <  —  2  *  10-17s-3m6-1.  Results  of  these 
two  sets  of  comparisons  are  shown  in  Tables  3  and  4.  Table  3  shows  the  ratios  while 
Table  4  shows  the  percentages. 

TABLE  3.  Ratio  of  number  of  precipitation  points  within  lifting  regions  to 
total  number  of  lifting  points  for  1200  UTC  17  March  1987.  Each  column  shows 
the  ratio  within  the  range  of  V  •  Q  given.  The  composite  column  considers  points 
with  lifting  at  any  level.  Units  of  V  •  Q  are  10-17s-3mfe 1 . 


800  mb  600  mb  400  mb  Composite 


382 

283 


TABLE  4.  Percentages  for  data  in  TABLE  3. 


mb  600  mb  400  mb  Composite 


400  mb 


The  results,  in  Tables  3  and  4,  show  that  using  V  •  Q  <  —2  gives  higher 
percentage  values  at  the  800  and  600  mb  levels  than  those  obtained  using  V  •  Q  <  0. 
In  this  case,  the  400  mb  level  did  not  show  a  significant  change.  The  composite 
results  were  determined  using  upward  vertical  motion  at  any  level  and  showed  a 
overall  increase  in  percentage  with  V  •  Q  <  —  2. 
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One  may  infer,  from  these  comparisons,  that  if  adequate  moisture  is  available, 
then,  areas  with  higher  values  of  upward  vertical  motion  have  an  increased  proba¬ 
bility  of  precipitation  occurring.  This  was  the  case  with  the  results  obtained  using 
values  of  implied  lifting  at  the  800  and  600  mb  levels. 

In  this  case  the  Barnes  program  showed  that  the  computed  quasi-geostrophic 
forcing  was  in  agreement  with  the  atmospheric  dynamics  at  the  time.  It  is  believed, 
based  upon  continuity  from  one  data  base  time  to  the  next,  that  the  NGM  was 
likely  to  continue  being  12  h  slow  in  its  precipitation  and  vertical  velocity  forecasts. 
However,  the  NGM  products  were  not  available  to  verify  this  supposition.  The 
Barnes  program  results  can  be  used  to  help  diagnose  the  current  state  of  the 
atmosphere  and  can  be  used  in  diagnosing  the  performance  of  other  products.  The 
Barnes  program  shows  more  detail  in  the  vertical  structure  of  the  atmosphere  than 
the  NGM  products.  Vertical  motion  is  implied  at  three  levels  with  the  Barnes 
program  while  the  NGM  only  gives  vertical  motion  at  one  level.  The  results  of  the 
Barnes  program  are  available  sooner  than  the  NGM  products.  The  Barnes  program 
results  are  available  approximately  10  min  after  acquisition  of  the  required  upper 
air  data. 


THUNDERSTORMS  IN  THE  CENTRAL  PLAINS  ON  4  JUNE  1980 


For  this  case  stud}'  the  Q-vector  charts  are  examined  for  correspondence  with 
a  low  level  jet  (LLJ)  and  the  initiation  of  a  large  scale  thunderstorm,  or  possible 
MCC,  event  over  central  Nebraska  and  also  over  northwestern  North  Dakota.  This 
is  done  in  order  to  check  the  Q-vector  results  under  a  different  synoptic  situation 
than  the  previous  case  study. 

The  synoptic  situation  at  0000  UTC  on  4  June  1980  had  a  quasi-stationary 
warm  front  extending  from  a  low  pressure  center  over  northwestern  Nebraska, 
southeastward  through  central  Missouri.  The  interest  in  this  case  study  is  in  the 
cold  air  ahead  of  the  warm  front  (Fig.  32). 

The  850  mb  geopotential  and  temperature  field  (not  shown)  showed  strong 
baroclinicity  with  the  isotherms  being  180°  out  of  phase  with  the  geopotential 
contours  throughout  the  central  plains  states.  The  wind  field  at  850  mb  showed 
a  LLJ  which  extended  from  the  panhandle  of  Texas  northward  through  central 
Nebraska  (shown  in  Fig.  32). 

The  700  mb,  500  mb,  and  300  mb  charts  (not  shown)  for  the  same  time 
showed  that  the  isotherms  at  these  levels  were  for  the  most  part  in  phase  with 
the  geopotential  contours.  Thus  the  atmosphere  was  equivalently  barotropic  at  and 
above  700  mb. 

The  thunderstorm  activity  over  North  Dakota  started  at  about  2000  UTC  on 
3  June  1980  while  the  activity  over  central  Nebraska  started  at  about  2300  UTC. 
These  times  were  determined  by  using  available  satellite  photographs  and  radar 
summaries.  By  0500  UTC  on  4  June  1980  each  of  the  storm  areas  covered  two 
states.  The  thunderstorms  over  North  Dakota  propagated  to  the  east-northeast 


FIG.  32.  850  mb  wind  field  and  surface  frontal  position  for  0000  UTC  4  June 
1980.  Wind  speed  is  in  knots.  Long  feathers  represent  10  kt  and  short  feathers 
represent  5  kt.  Surface  frontal  position  is  from  the  NMC  surface  chart  for  the  same 
time. 

while  the  Nebraska  storms  propagated  to  the  east. 

Referring  back  to  the  LLJ  at  850  mb  (Fig  32.),  the  maximum  wind  was  actually 
not  at  this  level.  Plotting  a  vertical  sounding  for  Oklahoma  City  showed  that  the 
maximum  winds  were  occurring  at  about  900  mb.  There  was  an  inversion  just  above 
the  850  mb  layer. 

It  is  seen  in  Fig.  33,  that  there  is  an  area  of  strong  Q-vector  convergence  near 
the  axis  of  the  southward  motion.  This  produces  the  area  of  maximum  upward 
motion  over  Nebraska  and  Kansas  seen  in  Fig.  34.  This  upward  motion  area 
coincides  with  the  Nebraska  thunderstorm  system.  Over  North  Dakota  there  is  a 
large  value  of  lifting  which  may  help  to  explain  the  development  of  the  thunderstorm 
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FIG.  33.  800  mb  Q  and  temperature  for  0000  FTC  4  June  1980.  Arrows 

represent  the  direction  and  magnitude  of  Q.  Temperature  is  given  in  °C. 
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FIG.  34.  800mb  V-Q  for  0000  UTC  4  June  1980.  Units  are  10  17  s  3  mb  1. 


The  600mb  V  Q  chart  (Fig.  35)  for  0000  I  !  (’  1  June  shows  that  there  was 


upward  motion  over  the  same  areas  but  that  the  magnitude  was  smaller.  From  this 
decrease  in  V  •  Q  over  North  Dakota,  Nebraska,  and  Kansas,  one  may  infer  that 
the  forcing  which  contributed  to  the  storm  environment  was  strongest  in  low  levels 
and  probably  associated  with  the  LLJ. 


FIG.  35.  600  mb  V  •  Q  for  0000  UTC  4  June  1980.  Labelling  is  the  same  as 

in  Fig.  34. 

The  800  mb  GFF  chart  (Fig.  36)  for  0000  UTC  4  June  shows  that  the  flow, 
depicted  in  Fig.  33,  has  enhanced  the  thermal  gradient  over  the  Minnesota  and 
Iowa  area.  The  strong  positive  GFF  area  is  on  the  cold  side  of  the  warm  front  and 
one  may  assume  that  the  circulation  in  that  area  was  favorable  for  maintenance  or 
enhancement  of  the  front. 


Twelve  hours  later,  at  1200  UTC,  both  storms  continued  to  grow7  in  size. 
The  800  mb  Q-vector  chart  (Fig.  37)  for  1200  UTC  4  June  1980  shows  that  the 
convergence  area  has  moved  eastward  over  the  preceeding  12  h.  The  convergence 
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FIG.  36.  800  mb  GFF  for  0000  UTC  4  June  1980.  Units  are  described  on  the 
bottom  of  the  map. 


area  is  now  over  the  western  Iowa  area  and  this  coincides  with  the  eastward 


translation  of  the  northern  portion  of  the  LLJ  during  the  same  period.  The  Q- 


vectors  pointing  southward  still  mark  the  area  of  the  LLJ. 


The  800  mb  V  •  Q  chart  for  1200  UTC  4  June  1980  (Fig.  38)  shows  the 
area  of  upward  motion  has  propagated  eastward  into  the  same  region  and  that  the 


magnitudes  of  V-Q  have  increased.  The  upward  motion  area  over  North  Dakota  has 


grown  in  size  while  the  subsidence  area  over  Minnesota  is  now  confined  to  northern 


Minnesota  and  has  weakened  over  the  12  h  period. 


The  800  mb  GFF  chart,  for  1200  UTC  4  June  1980,  (Fig.  39)  shows  that 


the  gradient  of  positive  values  over  Iowa,  eastern  Nebraska,  and  northern  Missouri 


has  increased.  From  this  increase  one  may  expect  that  the  front  now  has  a  more 


vigorous  circulation  and  is  more  intense. 


By  0000  UTC  on  5  June  1980  all  the  storms  discussed  were  diminishing  in 
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FIG.  39.  800  mb  GFF  for  1200  UTC  4  June  1980.  Labelling  is  the  same  as  in 
Fig.  36. 


size  and  intensity.  The  LLJ  at  this  time  was  also  weakening.  The  800  mb  Q- 
vector  chart,  for  this  time,  (Fig.  40)  shows  smaller  Q-vectors.  This  implies  that 
the  corresponding  flow  was  weaker  and  thus  the  vertical  motion  as  a  result  of  the 
convergence  of  Q  was  weaker.  With  a  decrease  in  lifting  the  storm  should  weaken. 

The  Q-vector  convergence,  shown  in  Fig.  41,  also  shows  smaller  values 
over  the  Nebraska  and  Iowa  area.  The  magnitude  of  the  convergence  increased 
over  Oklahoma  and  eastern  portion  of  the  Texas  panhandle.  Coincidently  a 
thunderstorm  developed  over  the  northern  Texas  area  and  lasted  for  a  few  hours. 

The  positive  GFF  (Fig.  42)  values  for  0000  UTC  on  5  June  1980  have  decreased 
in  magnitude.  This  was  primarily  due  to  the  smaller  magnitude  of  the  Q-vectors 
over  Iowa  (Fig.  40)  and  thus  there  was  weaker  ageostrophic  motion  over  Iowa. 

So  far,  this  case  study  has  shown  that  by  using  the  Q-vector  diaenostic  charts 
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FIG.  40.  800  mb  Q  and  temperature  for  0000  UTC  5  June  1980.  Labelling  is 
the  same  as  in  Fig.  33. 
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FIG.  41.  800  mb  V  •  Q  for  0000  UTC  5  June  1980.  Labelling  is  the  same  as 

in  Fig.  34. 


FIG.  42.  800  mb  GFF  for  0000  UTC  5  June  1980.  Labelling  is  the  same  as  in 
Fig.  36. 

an  operational  meteorologist  could  have  seen  the  changes  in  atmospheric  circulations 
and  made  an  improved  forecast  accordingly.  For  example,  the  increase  in  values  of 
the  upward  motion  and  also  GFF  could  have  led  the  forecaster  to  conclude  that  the 
thunderstorm  activity  would  have  continued  over  Iowa.  It  has  also  been  seen  that 
the  strongest  upward  vertical  motion  was  in  the  vicinity  of  the  LLJ. 

In  this  case  study  the  time  period  considered  was  from  0000  UTC  4  June  1980 
through  0000  UTC  5  June  1980.  In  order  to  check  the  correspondence  between 
different  combinations  of  precipitation  and  vertical  motion  for  both  the  mature  as 
well  as  the  decaying  stages  of  a  thunderstorm,  the  time  period  considered  will  be 
extended  to  0000  UTC  6  June  1980. 

A  set  of  composite  charts  was  constructed  using  the  same  methods  described 
in  the  previous  case  study.  The  resultant  composite  chart  with  V  •  Q  for  0000 


UTC  4  June  1980  and  the  NWS  Radar  Summary  for  0135  UTC  4  June  1980  are 
depicted  in  Fig.  43.  Figure  43  shows  that  we  do  not  have  a  large  continuous 
(homogeneous)  precipitation  area.  Here  we  have  relatively  small  and  somewhat 
separated  precipitation  areas.  This  is  also  true  for  the  composite  charts  created  for 
the  remaining  four  time  periods  under  study.  Each  time  period  is  12  h  later  that 
the  previous  one.  Since  these  charts  are  similar  to  Fig  43,  thay  are  not  shown.  We 
are  primarily  interested  in  the  ratios  evaluated  from  these  charts. 


FIG.  43.  Composite  chart  for  0000  UTC  4  June  1980.  Precipitation  (shaded 
area)  is  overlaid  with  computed  values  of  V  •  Q  <0  for  the  800  mb  (short  dashed 
lines),  600  mb  (long  dashed  lines),  and  400  mb  (dotted  lines)  levels.  The  arrows 
represent  the  side  of  the  contour  with  V  •  Q  <  0. 

Concerning  the  tables  that  follow,  the  Time/Date  column  shows  the  time  of 
the  Radar  Summary  chart  used.  The  Radar  Summary  charts  are  compared  with  the 
V  •  Q  charts  computed  for  the  nearest  time.  In  all  but  the  first  case,  in  each  table, 
there  is  a  25  min  difference  in  time  between  the  V  ■  Q  chart  and  the  corresponding 


Radar  Summary  chart.  In  the  first  case  there  is  a  95  min  time  difference.  The  2335 
UTC  3  June  1980  NWS  Radar  Summary  chart  was  not  available  for  some  unknown 
reason. 

The  first  set  of  comparisons  determines  the  ratio  of  number  of  precipitation 
points  with  implied  lifting  to  the  total  number  of  precipitation  points.  This  is  done 
for  the  800.  600,  and  the  40Q  mb  levels.  A  composite  ratio  is  also  determined  by 
considering  lifting  at  any  level. 

Table  5  shows  the  results  of  determining  these  ratios.  Table  6  shows  the 
percentages  for  the  data  in  Table  5.  The  800  mb  level  shows,  in  Table  6,  the 
highest  percentages,  of  precipitation  points  having  implied  lifting  motion,  in  four 
out  of  five  cases.  The  600  mb  level,  in  three  out  of  five  cases,  shows  a  relatively 
large  drop  in  percentages  from  values  obtained  at  the  800  mb  level.  In  all  five 
cases,  the  400  mb  level  had  the  lowest  percentages.  For  the  composite  column,  the 
percentages  remained  fairly  high.  At  this  stage  we  may  infer  that  the  precipitation 
was  being  caused  by  primarily  lower  to  middle  tropospheric  effects. 

TABLE  5.  Ratio  of  precipitation  points  with  lifting  over  total  number  of 
precipitation  points  for  case  2.  The  composite  column  considers  points  that  have 
lifting  at  any  level. 


Time/Date 

800  mb 

600  mb 

400  mb 

Composite 

0135  UTC  4  JUNE 

29/41 

30/41 

20/41 

35/41 

1135  UTC  4  JUNE 

33/50 

22/50 

21/50 

37/50 

2335  UTC  4  JUNE 

52/72 

32/72 

25/72 

60/72 

1135  UTC  5  JUNE 

32/51 

30/51 

20/51 

36/51 

2335  UTC  5  JUNE 

33/55 

24/55 

18/55 

46/55 

For  this  case  study,  we  are  not  se^  ng  the  same  consistency,  between  levels, 
that  we  saw  with  the  previous  case  study.  Part  of  the  reason  may  be  due  to  the 


ivivj  vr  *»;« 


71 

73 

66 

44 

72 

44 

63 

59 

44 

TABLE  6.  Percentages  for  data  in  TABLE  5. 


Time/Date  800  mb  600  mb  400  mb  Composite 


0135  UTC  4  JUNE 
1135  UTC  4  JUNE 
2335  UTC  4  JUNE 
1135  UTC  5  JUNE 
2335  UTC  5  JUNE 


geographical  location  of  the  precipitation  in  the  northwestern  portion  of  the  grid. 
Notice,  in  Fig.  43,  that  the  northwestern  portion  of  the  grid  is  located  over  Idaho 
and  western  Montana.  This  area  contains  a  portion  of  the  Rocky  Mountain  range. 
Therefore  some  of  the  precipitation,  in  this  area,  may  be  due  to  terrain  effects  and 
not  due  to  quasi-geostrophic  effects.  Referring  back  to  the  composite  chart  (Fig. 
43),  it  can  be  seen  that  if  the  precipitation  areas  over  the  northwestern  area  of 
the  grid  are  disregarded,  then  the  percentage  of  precipitation  areas  having  upward 
vertical  motion  would  increase  dramatically  at  every  level.  Since  it  is  known  that  the 
atmosphere  does  not  react  only  to  quasi-geostrophic  effects,  the  precipitation  area 
in  question  will  not  be  disregarded.  In  all  five  time  periods  discussed,  the  majority 
of  the  precipitation  occurring  within  subsidence  areas  is  in  the  Idaho/Montana  area. 

The  confidence  interval  (99%)  for  a  binomial  experiment  was  evaluated  for  the 
800  mb  values  shown  in  Table  6.  In  every  case  the  percentages  shown  in  Table  6 
were  found  to  be  statistically  significant.  This  was  expected  since  precipitation  is 
usually  associated  with  low  level  upward  motion. 

Next  the  ratio  of  the  subsidence  areas  that  are  precipitation  free  is  examined. 
For  this  comparison,  Table  7  shows  the  tabulated  data  while  Table  8  shows  the 
percentages  for  the  data  in  Table  7.  The  composite  column,  in  each  of  these  tables, 
the  subsidence  is  determined  in  areas  that  show  subsidence  at  every  level. 
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K  Time/Date 

800  mb 

600  mb 

400  mb 

Composite 

P  0135  UTC  4  JUNE 

95 

96 

93 

97 

L 

!  1135  UTC  4  JUNE 

92 

86 

88 

89 

;■  2335  UTC  4  JUNE 

92 

87 

85 

93 

?  1135  UTC  5  JUNE 

93 

92 

89 

92 

V 

I  2335  UTC  5  JUNE 

92 

97 

95 

95 

TABLE  7.  Ratio  of  number  of  subsidence  points  that  are  precipitation  free  over 
the  total  number  of  subsidence  points  for  case  2.  The  composite  column  considers 
points  with  subsidence  occuring  at  every  level. 


Time/Date 

800  mb 

600  mb 

400  mb 

Composite 

0135  UTC  4  JUNE 

241/253 

259/270 

287/308 

177/183 

1135  UTC  4  JUNE 

205/222 

170/198 

216/245 

108/121 

2335  UTC  4  JUNE 

225/245 

268/308 

259/306 

157/169 

1135  UTC  5  JUNE 

268/287 

256/277 

263/294 

180/195 

2335  UTC  5  JUNE 

223/242 

270/279 

322/340 

155/164 

TABLE  8.  Percentages  for  data  in  TABLE  7. 


Table  8  shows  that  the  results  remained  relatively  constant  through  each  time 
period  and  each  level.  These  results  are  very  similar  to  the  corresponding  results  in 
the  previous  case  study  (see  Tables  1  and  2).  It  can  be  inferred  that  the  computed 
subsidence  areas  are  relatively  accurate  since  these  areas  are  almost  precipitation 
free,  as  expected. 

Two  more  sets  of  comparisons  are  made  to  conclude  this  study.  Both  sets 
determine  the  ratio  of  precipitation  points  within  a  area  of  lifting  to  the  total 
number  of  implied  lifting  points.  The  difference,  in  these  two  sets,  is  the  value  of 
V  -Q  used.  The  first  set  considers  vertical  motion  with  computed  values  of  V-Q  <  0 
while  the  second  set  considers  values  of  V  ■  Q  <  —2  *  10_1 '  s~3mb~^ . 

The  results  using  V  •  Q  <  0  are  shown  in  Tables  9  and  10.  Table  9  shows  the 
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evaluated  data  and  Table  10  shows  the  percentages  of  the  data  in  Table  9.  For  these 
tables,  the  composite  column  considers  lifting  at  any  level.  As  seen  in  these  tables 
the  ratio  of  the  total  lifting  area  having  precipitation  is  relatively  small. 

TABLE  9.  V  •  Q  <  0  versus  precipitation.  Each  column  shows  the  ratio  of  the 
number  of  lifting  points  having  precipitation  to  the  total  number  of  lifting  points. 
The  composite  column  considers  lifting  at  any  level. 


Time/Date 

800  mb 

600  mb 

400  mb 

Composite 

0135  UTC  4  JUNE 

29/277 

30/260 

20/222 

35/347 

1135  UTC  4  JUNE 

33/308 

22/332 

21/285 

37/409 

2335  UTC  4  JUNE 

52/285 

32/222 

25/224 

60/361 

1135  UTC  5  JUNE 

32/243 

30/253 

20/236 

36/335 

2335  UTC  5  JUNE 

33/288 

24/251 

18/190 

46/366 

TABLE  10.  Percentages  for  data  in  TABLE  9. 


Time/Date 

800  mb 

600  mb 

400  mb 

Composite 

0135  UTC  4  JUNE 

10 

12 

9 

10 

1135  UTC  4  JUNE 

11 

** 

i 

7 

9 

2335  UTC  4  JUNE 

18 

14 

11 

17 

1135  UTC  5  JUNE 

13 

12 

8 

11 

2335  UTC  5  JUNE 

11 

10 

9 

13 

Tables  11  and  12  show  the  ratio  of  the  total  lifting  area  having  precipitation. 
Lifting  values  used  for  this  comparison  are  V  •  Q  <  —2  *  10_17s_3m6_1 .  Table  11 
shows  the  evaluated  data  wrhile  Table  12  shows  the  percentages  for  the  data  in  Table 
11.  Comparing  Table  10  with  Table  12  shows  that  the  800  mb  percentages  increased 
with  V  •  Q  <  —2  *  10  11  s  3 mb  1 .  This  is  also  true,  in  three  out  of  five  cases,  for 
the  600  mb  level.  The  400  mb  level  showed  a  drop  in  percentages  with  the  stronger 
lifting  value  in  all  five  cases.  Finally,  the  composite  column  shows  a  increase  in 
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TABLE  11.  V-Q<—  2*10  1 ' s  3mb  1  versus  precipitation.  This  table  is  the 
same  as  Table  9  except  for  the  value  of  V  •  Q  used  in  determining  the  ratios. 


Time/Date 

800  mb 

600  mb 

400  mb 

Composite 

0135  UTC  4  JUNE 

18/137 

12/106 

7/112 

22/204 

1135  UTC  4  JUNE 

17/154 

7/135 

9/192 

20/290 

2335  UTC  4  JUNE 

24/104 

10/55 

4/51 

28/156 

1135  UTC  5  JUNE 

30/87 

21/77 

0/77 

30/152 

2335  UTC  5  JUNE 

28/165 

14/104 

11/76 

29/199 

TABLE  12.  Percentages  for  data  in  TABLE  11. 


Time/Date 

800  mb 

600  mb 

400  mb 

Composite 

0135  UTC  4  JUNE 

13 

11 

6 

11 

1135  UTC  4  JUNE 

11 

5 

5 

7 

2335  UTC  4  JUNE 

21 

18 

8 

18 

1135  UTC  5  JUNE 

34 

27 

0 

20 

2335  UTC  5  JUNE 

17 

13 

6 

15 

percentages  for  four  out  of  five  cases.  It  is  again  inferred  that  the  precipitation  may 
have  been  due  to  primarily  lower  to  middle  tropospheric  forcing. 

One  can  also  imply,  from  Tables  10  and  12,  that  given  adequate  water 
vapor,  quasi-geostrophic  effects  in  the  lower  troposphere  are  more  likely  to  produce 
precipitation  than  the  same  effects  in  the  upper  troposphere.  Normally,  an 
operational  meteorologist  will  look  for  adequate  water  vapor  and  lifting  to  produce 
precipitation.  The  meteorologist  usually  looks  for  these  parameters  in  the  lower  to 
middle  troposphere.  Therefore,  it  is  seen  that  the  above  results  are  consistent  with 
concepts  commonly  used  by  a  meteorologist. 
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CHAPTER  V 


CONCLUSIONS 


This  paper  gives  a  review  of  the  Barnes’  diagnostic  model  and  its  application  in 
operational  work.  The  examples  presented  demonstrate  the  usefulness  of  Q,  V  •  Q, 
and  GFF  in  a  variety  of  weather  situations. 

In  the  case  of  the  thunderstorms  through  eastern  Texas,  the  Barnes’  diagnostic 
charts  outperformed  the  much  more  elaborate  NGM  products.  The  Q-vector,  V Q, 
and  GFF  charts  correctly  showed  the  quasi-geostrophic  forcing  for  the  frontal  system 
moving  through  Texas.  The  NGM  products  were  approximately  12  h  slow  in 
movement  of  the  front  and  associated  precipitation.  The  forcing  depicted  in  the 
Barnes  program  results  was  also  conducive  to  enhancement  of  the  thunderstorms 
and  subsequent  severe  weather.  This  case  study  also  showed  additional  charts 
produced  by  the  Barnes  program  that  may  be  of  use  in  operational  meteorology. 
These  include  mixing  ratio,  geopotential  tendency,  and  stability  tendency  charts. 
Mixing  ratio  charts  are  computed  for  four  levels  in  the  atmosphere  while  the 
geopotential  and  stability  tendency  charts  are  computed  for  three  levels.  This  gives 
the  forecaster  a  more  detailed  picture  of  the  vertical  structure  of  the  atmosphere. 

The  second  case  study  showed  the  location  and  subsequent  spreading  of  an 
MCC  in  detail.  The  depicted  forcing  agreed  with  thunderstorm  and  front  locations. 
In  this  case  the  low  level  forcing  was  shown  to  be  more  significant  than  middle 
and  upper  level  forcing.  The  800  mb  V  •  Q  chart  showed  the  strongest  convergence 


magnitudes  and  these  occurred  in  the  vicinity  of  the  LLJ. 

For  both  case  studies  the  precipitation  field  was  for  the  most  part  within  an  area 


of  computed  lifting.  The  computed  subsidence  areas  were  primarily  precipitation 
free.  This  simple  check  lends  credence  to  the  computed  V  •  Q  values.  Also  shown 
in  both  cases  was  that  precipitation  was  more  likely  to  occur  as  the  magnitude  of 
the  convergence  of  Q  increased  in  the  800  and  600  mb  levels.  GFF  gave  a  good 
indication  of  frontal  intensity  and  subsequent  development. 

Overall,  the  Barnes  program  provides  a  more  detailed  picture  of  the  quasi- 
geostrophic  effects  in  the  atmosphere.  The  Q-vector  charts  give  information  on 
the  horizontal  ageostrophic  motion  which  in  turn  gives  an  indication  of  the  frontal 
circulations.  The  Q-vector  field  overlaid  with  a  layer  temperature  field  graphically 
shows  where  the  thermal  field  is  being  enhanced  or  weakened  and  thus  gives 
a  reliable  picture  of  frontal  development.  The  magnitude  of  the  Q-vector  is 
proportional  to  the  strength  of  the  horizontal  ageostrophic  wind  below  the  level 
at  which  Q  is  computed. 

The  Q-vector  formulation  of  the  omega  equation  and  the  corresponding  charts 
are  valuable  tools  when  used  in  conjunction  with  other  diagnostic  and  prognostic 
products.  As  is  well  known,  there  can  be  differences  between  different  charts.  Using 
a  program  such  as  the  Barnes  program  gives  the  meteorologist  another  tool  to 
help  ascertain  the  current  state  of  the  atmosphere.  Once  the  current  state  of  the 
atmosphere  has  been  determined  then  the  product  differences  can  be  resolved  and 
a  more  accurate  forecast  can  be  made. 

Two  potential  problems  arise  when  comparing  the  NWS  Radar  Summary  charts 
with  the  V  •  Q  charts.  First  the  radar  summary  charts  differ  in  time  from  the  V  •  Q 
charts  by  about  30  min.  Second,  the  radar  summary  charts  are  an  estimate  of  the 
actual  precipitation  fields.  The  radar  summary  charts  are  a  composite  of  the  many 
radar  reports  received  by  the  NWS  and  some  of  the  lighter  precipitation  falling  may 
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not  be  seen  by  the  radar. 


Future  research  efforts  could  possibly  use  the  actual  radar  reports  to  plot  the 
areas  of  precipitation  and  then  make  the  same  types  of  comparisons.  Another 
research  recommendation  is  to  determine  the  relationship  between  the  amount  of 
atmospheric  water  vapor,  V-Q  values,  and  precipitation.  The  final  recommendation 
is  for  a  frequency  distribution  study  of  the  values  of  V  •  Q  and  precipitation.  This 
frequency  distribution  study  may  provide  a  critical  value  of  V -Q  that  infers  probable 
precipitation.  This  would  be  similar  to  the  relationship  between  stability  indices 
such  as  K -Index  or  Total-Totals  values  and  convective  activity. 
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